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A standard perturbative expansion around the mean-field solution is used to derive the low-energy 
effective action for superconductors at T — 0. Taking into account the density fluctuations at the 
outset we get the effective action where the density p is the conjugated momentum to the phase 9 of 
the order parameter. In the hydrodynamic regime, the dynamics of the superconductor is described 

CI 

by a time dependent non-linear Schrodinger equation (TDNLS) for the field ^(x) = \f\e %e . The 
evolution of the density fluctuations in the crossover from weak-coupling (BCS) to strong-coupling 
(Bose condensation of localized pairs) superconductivity is discussed for the attractive Hubbard 
model. In the bosonic limit, the TDNLS equation reduces to the the Gross-Pitaevskii equation for 
O ' the order parameter, as in the standard description of superfluidity. The conditions under which a 

»' ! , phase-only action can be derived in the presence of a long-range interaction to describe the physics 

O |i of the superconductivity of "bad metals" are discussed. 
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INTRODUCTION 

■ 

G ■ 

The description of a superconducting system at low temperatures by an effective action acquired a renewed interest 
O | due to the recent attention to vortex dynamicsEJ and to the non-standard behavioiH of the High T c Superconductors 
(HTCS) . The advantage in having an effective theory is well known: the complex microscopic theory is replaced by a 
y—( | simpler one, written in,tcrms of relevant variables only. As an example, near the critical temperature, the equilibriumtl 
^ i and non-equilibriumQ#Q properties of a superconductor are well described by the Landau- Ginzburg theory that is 
I/"} ' written in terms of the complex order parameter A(x) = (cpi(x)(pr(x)). 

At zero temperature, the effective theory which allows a minimal description of the low-energy properties of a 
superconductor (from the Meissner effect to the flux quantization) is written in terms of the phase 9{x) of the order 
parameter. Microscopic derivations of such an effective theory are usually carried out by deriving the low temperature 
effective action in terms of the complex order parameter. This effective action then leads to a wave-egusution for the 
\ phase while the order parameter itself obeys a non-linear equation with second order time derivativesB'cra. r— ■ 
By contrast, the low-energy dynamics of a bosonic superfluid is well described by a Gross-Pitaevskii equationt2l 
i.e. a time-dependent non-linear Schrodinger equation (TDNLS) for the complex order parameter \P = |^/;|e l6 '. The 
continuity equation implied for |^(a;)| 2 by the TDNLS description is a continuity equation for the full density which, 
in this case, can indeed be identified with |^(a;)| 2 . For a fermionic superfluid there is instead a non-trivial relation 
between the density p(x) and the order parameter A(x). However, the Gross-Pitaevskii description is expected to 
hold also for fermions in the strong pairing limit, when the Cooper pairs, that are strongly overlapping in the BCS 
\ description, can be considered as composite bosons. In this limit, the phenomenological effective lagrangian proposed 
• • ■ by Feynmantj, suggesting description o£-the superconducting system as a bosonic superfluid, should be recovered. 

Recently, Aitchson et al&3 and M.StoneEj re-interpreted the wave-like equation for the phase of the order parameter 
for a BCS superconductor in terms of a non-linear Schrodinger equation. They reconciled the wave-equation with a 
TDNLS description by re- writing the low-energy effective lagrangian for A(x) in terms of the phase and of the field that 
plays the role of the conjugated momentum of 9(x) i.e. the local density p(x). The field p(x) is however introduced 
in the theory as an ad hoc field. This procedure needs indeed an identification a posteriori of the "wave-function" 
obeying the non-linear Schrodinger equation and does not explain the strong coupling limit. 

Here we give a microscopic derivation of the zero-temperature low-energy effective action in terms of relevant 
physical fields, including from the beginning the density p(x) together with the complex order parameter A(x)B. Our 
derivation is quite general. We consider, as usual, small fluctuations of the fields from their mean field values and 
retain in the effective action only the leading momentum and frequency dependences (hydrodynamic limit). 

The effective action is obtained by a standard perturbative procedure!!! with two main differences with respect to 
previous approaches. First, we decouple the pairing term by means of an Hubbard-Stratonovich (HS) transformation 



X 



1 



in the particle-particle and in the particle-hole channels to include both fluctuations. Since the auxiliary fields 
introduced by the Hubbard-Stratonovich transformation are not the physical fields, the latter fields are introduced 
by a further transformation. This makes the low-energy effective lagrangian easily readable and allows for a direct 
interpretation of the local density as the conjugated momentum of the phase. Furthermore the correct strong pairing 
limit is straightforward. 

The equations of motion supported by this microscopic derivation are the non-linear Schrodingcr equation for the 
variable ^(x) = \J ^^-e t9( - x \ in agreement with Refill andEl, and an equation for the modulus of the order parameter 
A(x)| relating the local condensate to the local density. 

When we consider the BCS-Bose crossover for the attractive Hubbard model, the expression of the lagrangian in 
terms of the physical fields shows that, on increasing the pairing interaction, the density tends to fluctuate coherently 
with |A(x)| 2 . In particular, in the low-density limit, the two fields become proportional and the description of low- 
energy dynamics by means of two equations of motion becomes redundant. The low-energy dynamics is then described 
by a Gross-Pitaevskii equation (a TDNLS equation for the complex order parameter A (a;)). 

On various occasions the crossover between the BCS (weak-coupling) and the Bose (strong-coupling) limit has been 
studiedEjES and the analysis of this problem has acquired particular relevance for the High T c Superconductors. 
Indeed, some common characteristics-of these compounds (the small value of fc^C — 10 — 20, compared to the BCS 
values, and the small carrier densityEj) can be in favor of a situation close to condensation of localized Cooper pairs, 
especially in the under-doped HTCS. In this context, the role of the phase fluctuations for -this compounds has been 
analyzed looking at the phase-only action in the presence of a long-rang e Coulomb forceMEi The low -energy action 
which we obtain, including the density fluctuations together with the fluctuations of the complex order parameter, 
easily allows to study situations that can deviate from the BCS behavior. Therefore, at the end, we also considp; 
long-range forces in our approach and compare (where possible) our result with the Emery and Kivelson proposaO 
of a phase-only effective action for a description of HTCS as "bad metals". 

The paper is organized as follows. The first section is dedicated to the derivation of the low-energy effective 
lagrangian in terms of the physical fields for a simple continuum model. In section II, we derive the equations of 

motion and obtain the non- linear Schrodinger equation for the field \P(x) = 



In the third section we 



consider the explicit crossover from weak to strong coupling starting from the attractive Hubbard model. Finally, the 
effects of a long-range interaction are specifically considered in the last section where the phase-only action is also 
derived. 

I. THE EFFECTIVE ACTION 

The details of the interaction inducing pairing are not relevant for our purposes. We then consider a very simple 
model: a fermionic gas, in 3-d dimensions, with a local pairing interaction and an energy cut-off to remove ultraviolet 
divergences. The "partition function" at T = (h = ks = 1) for this system is: 



V(j) + (x)V(j)(x)e tSi4 ' +ix) '^ x}) 
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(1.1) 



where x = (t,r) ,<^>+(x) and <p a (x) are the Grassmann variables associated to fermions of spin a, ji is the chemical 
potential and U < is the attractive interactions. 

To obtain the effective action in terms of relevant variables for the low-energy region, we decouple the interaction 
term by means of an HS transformation in both the particle-particle and the particlc-hole channel, introducing in this 
way the auxiliary fields associated to the complex order parameter and to the density. The decoupling of the pairing 
term reads: 
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prescription: 
The relative 



particle-hole channel 

weights a and 7 of the p-p and p-h channels are arbitrary but. for the "completeness" condition a 2 +7 2 =l. All the 
resulting effective actions are equivalent when treated exactlyEj allowing that each channel contributes to the other 
one at large q. However,_tp have the mean-field results as a starting point, we disregard the "completeness" condi- 
tion and take a = 7 = lE3 while introducing a cut-off in the momentum space to prevent over-counting. Indeed, a 
restricted integration on q < q* separates the contributions related to the two channels. Moreover, from now on, we 
will not consider the field p sp i n (x) since, in the low-energy limit, spin density fluctuations are completely decoupled 
from the rest of the system. 

In the following we briefly report the main steps of the derivation of the effective action, while its full derivation 
is reported in the Appendix. As a first step, since the HS felds p HS , A" s , (A HS ) + act as conjugate to the physical 
fields and coincide with them only at the saddle point levelEj, we introduce the physical fields by means of a further 
"Gaussian" transformatimi of the quadratic terms appearing in (1.2). Secondly, to explicit the dependence on the 
gradients of the phase 6@ll3 of order parameter in the action ( Al), we make a gauge transformation 4>(x] — » cj)(x)e lS ( x '' 2 . 

To obtain the effective action we retain the first and the second order of a standard expansion!] of the effective 
bosonic action (A3) which is obtained integrating out the fermionic fields, around the saddle-point solution in the 
superconducting state. As usual, we assume that the fluctuations of the fields around their saddle point values are 
small and in particul ar th at the variations of the phase in space and time are slow. The last step consists in integrating 
out the HS fields in (|A6|) . 

The effective action in terms of the physical fields reads: 
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The explicit expressions of the above coefficients at T=0 are reported in the Appendix. 

It is easily checked that eq.(1.3) reproduces the generalized RPA approximation for the response function in the low- 
frequency and low-momentum limitt3. The low-energy and low-momentum excitations of a superconducting sy stem 
described by a microscopic action like (1.1) are density fluctuations: i.e. the Anderson-Bogolubov sound modeEBEia. 
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I n or der to describe this collective mode, it is sufficient to consider only the hydrodynamic limit of the effective action 
(1.3) by taking the q, ui — > limit of the coefficients of the quadratic part in p and |A|. No tice that at zero temperature 
the static (us = 0, q— ► 0) and the dynamic limit (q = 0, us — ► 0) of the coefficients in fll.3|) coincides. 



In the BCS limit we recover the usual dispersion for Anderson-Bogolubov mode 



[1 + UN ] q 2 /3 = 



where vf — kp/m is the Fermi velocity and Nq is the density of states per spin at the Fermi energy. Notice that, 
in the BCS limit, 8p(x) and <J|A|(x) are decoupled since 11% = 0. Indeed, Bo oc J deN(e)e vanishes due to the 
particle-hole symmetry property of the BCS (weak-coupli ng) theory. Contributions beyond the hydrodynamic limit 
(i.e. a gradient expansion of the coefficients in the action ( |l.3j )) represent higher order corrections to the sound mode 
dispersion. 

Making a direct comparison with previous derivationsliirlij we find that the density po in front of the phase term has 



d t e(x)/2 
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are eliminated 



become the full local density p(x) — 5p(x) + po? while higher order contributions in 
away under the integration of the HS fields. A direct way to obtain this result is to make an appropriate translation 
of fields in the action ( |A3| ) which is expressed in terms of both the HS fields and the physical fields. We use the 
following transformation of unitary Jacobian: 



p HS (x) = p HS (x) + [d t 6(x)/2 + (\7d(x)) 2 /8m] 
and rewrite the term p HS (x)p(x) in the action ( |A3| ) as p HS (x)p(x) = {p HS (x) 



(1.4) 

(xj6(x)Y /Sm}} p(x). The 



second term of this expression is recognized as the first contribution in the action fll.q ) . Then the integration over 
p HS (x) and the other HS fields (within the same approximations used previously) results into the phase independent 
contribution of ( jl.3| ). By means of this procedure, we see that the first term of (1.3) is always present in the full 
action. This result strongly relies on the introduction of the density fluctuations at the outset of our approach and 
on the use of the physi cal fi elds p and A. 

The transformation (1.4) also allows us to see that the dependence of the action on dt0(x)/2 shown in eq.( |l.3| ) is 
not restricted to the Gaussian approximation for the HS fields showing that p(x) is the proper coefficient of dt0(x)/2. 
This approach has to be contrasted to the one of Ref.(ll), where the density p(x) is deriveiLphenomenologically by 
reconstructing p(x) = po + Sp(x) via the Galilean invariance from mean field like equations.LJ 



II. THE LOW-ENERGY EFFECTIVE ACTION: THE NON-LINEAR SCHRODINGER EQUATION 



Let us consider, for completeness, the system in the presence of an electro-magnetic field A = (ip(x),A(x)) that 

d t 9/2 + eip, yO/2 -> \/9/2 - -A and 



we introduce in (1.3) through the minimal substitution on the phase dtO/2 
adding the free electromagnetic action. The effective action then reads: 
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V 9(x) 



[U° pp Sp(x) 2 + H° AA S\A(x)\ 2 + 26p(x)n° Ap 6\A(x)\ 
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(W(z) + -d t A(x)f + ( V A A(x)f 



(2.1) 



where we call n^, H° pA , T1 AA the q — > limit of the coefficients previously defined. 

We have already discussed the origin of t he fi rst term in ( |l.3| ). Because of this term the role of physical density is 
directly stated: the functional derivative of (2T) with respect to d t 9(x) gives —p(x)/2 and identify half of the physical 
density as the conjugated momentum to the phase. 

This feature suggests to introduce the field ^(x) — J £^£l e l0 ( x ) and to re-write an effective action in terms of it. 



The eq. (2T) becomes S e ff=J d 4 xC e ff with: 
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( V ^+-a f A) 2 + ( V AA)' 
c 



(2.2) 



Notice that the presen ce o f a ionic background, as usually assumed to conserve the charge neutrality of system, leads 
to the replacement in (2.2) of 2e(p(x)^*^ in e(p(x) (2\I r *\f r — pi on )- This term vanishes on average since pp — pi on , but 
it becomes eip(x) (2W*^ — po) by retaining the dens ity fl uctuations. 

The mapping of the low-energy effective action (2.1) into (2.2) holds in the strict hydrodynamic limit. Indeed 



only the term ^- (\/0) contained in y ^* y \& is present in the lagrangian (2.1), while the remaining term 



32mp (yp) should be evaluated from the gradient expansion of the quadratic contribution oc (Sp) in (1.3). However 
the coefficient of (yp) 2 we obtain from (1.3) is not 32 ] np (yp) 2 (for instance, this evaluation at half-filling, for the 
Hubbard model, gives the same expression §^-g jq we would obtain from expanding the compressibility at low q in 
the normal phase). Therefore the mapping is reliable only when this difference is irrelevant: i.e. in the low- frequency 
and low-momentum limit for the the density fluctuations. 

The introduction of the field ^(x) allow us to derive the equation of motion for ^(x) in form of a non-linear 



Schrodinger equation. By taking the functional derivatives of (2.2) with respect to ^>*, <5|A| and A we get: 

idtV(x) - 2etp(x)y(x) --?-( —A(x) + isj] #(x) - -(V(x)**(x) - p°/2)*(x) = 

4m V c / x 
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where x represents the compressibility given by: 
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(A - Dp) [1 + U (Dp + Cp)] - WBl 
(n° A ) 2 "[l + U{Ao- Dp)] [l + U (Dp + Cp)} - AU 2 Bp 



= 2 



(2.4) 
(2.5) 



and the current j(x) obeys the Maxwell equation: y A B = —j(x) + \dtE. 

Going back to the variables p(x) and 9(x), the imaginary and real part of the non-linear Schrodinger equation 



(2.3) are the continuity equation and the equation that links fluctuations in the density to fluctuations in the phase 
respectively: 



d t p(x) + y 



m \ 2 c w 







(2.6) 



Sp(x) = -x < [d t 9(x)/2 - e<p(x)] + — 



(2.7) 



The continuity equation (p.6|) is equivalently obtained by setting to zero the functional derivative of (2.1) with respect 



to 9(x). The appearance of the continuity equation is indeed due to the fact that the phase is the Goldstone boson 
field whose associated equation of motion is the charge conservation, provided the electromagnetic response is only 
given in terms of the phase (as it is the zero temperature case). 

As already pointed out, the difference (in general) between the modulus squared |A| 2 of the order parameter and the 
density p in a superconductor prevents to obtain a descript ion of the dynamics uniquely in terms of a Gross-Pitaevskii 
equation. In fact, the non-linear Schrodinger equation (|2.3| ) is related to the continuity equation and cannot be 
extended to describe the condensate |A|, but for the bosonic limit which we discuss in the next section. To obtai n th e 
dynamics of the condensate we should have made a gradient expansion of the coefficients of the effective action ( |l.3| ) , 
at least up to the second derivatives (in time and space) of the modulus |A| to define the coherence length for this 
field. Being our description limited to the hydrodynamic low-energy regime, the dynamics of the amplitude of the 
order parameter is missing. Indeed, the local eq.(^4|) should be more properly understood as a static relation, valid 
at q = 0, that links the density fluctuations to the order parameter amplitude fluctuations. In conclusion for almost 
constant |A|, the relevant dynamics is described by the non- linear Schrodinger equation for ^>(x). The coherent length 
we get for \l/(x) is linked to the compressibility and is then roughly of the order of inter-particle spacing. 
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III. THE BOSONIC LIMIT: THE ATTRACTIVE HUBBARD MODEL 



We want to analyze the problem of the crossover from standard BCS (weak-coupling) superconductivity to strong 
coupling superconductivity, where the Cooper pairs are made of couples of strongly bound fermions, by extending the 
previous approach to the attractive Hubbard model on a latticeEj. Our aim is to understand the interplay between 
density and order parameter while increasing the strength of the pairing interaction. 

In particular, we want to see how the description of the dynamics of a strong coupling superconductor evolves from 
the BCS limit to the Gross-Pitaevskii equation and how the role played by the local density is taken by the square of 
the order parameter amplitude, at least in the low-density limit. 

The fermionic action for the Hubbard model is: 

S(4> + (x),c/>(x)) = (idt + H) 4>i,a +tJ2 ««A* + #<A<0 - U Y, <T&.T<^U (3- 1 ) 

i,a (ij) i 

where, as in the previous section, U < 0, the sum (ij) is over the nearest neighbour pairs and the lattice spacing is 
a = 1. 

In deriving the effective action for the Hubbard model on a lattice two cases are of relevance: 

i) the hydro dynamic limit where one can assume small variations \0i — 6j\ ~ a\ V 9i\ <C 1, for the phase of the 
order parameter. In this case we follow the same procedure applied in the continuum limit, in particular we expand 



-iTrlnG 1 in (A3) in powers of 1 — cos — Oj) — — Oj) 2 to get the gradient contribution (v^i) 2 



ii) away from the hydrodynamic limit to allow \9i — 0j \ — 1. In this case the differences in the nearest neighbor site 
phases are not small and the expansion parameter \6i — 9j\ should be substituted by t/\U\. Indeed, it is only in the 
limit t/U <C 1 that one can assign a consistent meaning to the site phas e 9{l3. We shall consider this case in the next 



section, in connection with the derivation of the phase-only action (4.1). 

Here, to analyze the bosonic limit, we hold on the more standard hydrodynamic approach. The derivation is 
straightforward and the effective low-energy and low-momentum action in terms of the physical fields reads: 



S eff 



dx A [p(x)d t e(x)/2 - 
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[5 P (x)n° pp Sp(x) + S\A(x)\U° AA S\A(x)\ + 26 P (x)U Ap 6\A(x)\] } 



(3.2) 



n° pA ,H AA are the lattice version of the previously defined coefficients. Obviously, apart from differences 



where 11° , 

due to the lack of Galilean invariance, the effective lagrangians in the continuum (1.3) and in the lattice case ( |3.2| ) 
are similar. Indeed, as far as the explicit time dependence of action is concerned, we obtain the same combination 
p(x)dt9(x)/2 in both cases. The stiffness is here given in terms of the averaged kinetic energy (E^ m ), consistently 
with the expression of the diamagnetic contribution for the current, which changes from —po-§^ into §^(-E<f m )- 

In the attractive Hubbard model, except at half-filling where particle-hole symmetry holds, the coefficient, propor- 
tional to Bq, that links fluctuations of the modulus of the order parameter to the density fluctuations is finite and its 
role is essential for the description of the system in the strong pairing limit. 

In order to get informations on the asymptotic (large \U\) behaviour, we diagonalize the part of the lagrangian 
(3.2) which is quadratic in the fields Sp and S\A\. We obtain, in this way, two eigenvectors: 
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with eigenvalues: 
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(n°, p -A0 



A 1)2 = ~ {(Tl pp + H AA ) T [(H° pp n °AA) 2 + 4(H° A ) 2 ] V2 } 



On increasing the interaction, one of the eigenvalues (A2) diverges like U 2 , while the value of the other (Ai) stays 
finite. This behaviour implies that the eigenvector corresponding to the diverging eigenvalue, V2, has vanishing 
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fluctuations [v% V2) oc A 2 . Therefore Sp tends to fluctuate coherently with S\A\ through the coefficients of the linear 
combination appearing in V2'- 



S\A\ 



n A z 



As 



n 



pA 



Sp. 



(3.3) 



with the ratio ^ 2 staying finite for large U, since both IIaa and n p A diverge like U 2 . The density and the 
modulus of the order parameter experience the same fluctuations and we expect that they play the same role, so that 
a description of the low-energy dynamics of system by means of two equations of motion for Sp and S\ A| is becoming 
redundant. To clarify further the evolution of fluctuations of | A| we integrate out the density fluctuations and replace 
2|A |<5|A ~ <5|A| 2 . The resulting effective action is: 



S(0,S\A\ 



dx 



I TpKvn\ 
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MW) 2 -/ dx 4 {^d t 6(x) + 
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2A n 



ir 



AA 



11 Ap 11 pA 

n° 

pp 
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2A n 
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On increasing the interaction, independently from the filling, the term proportional to dtO(x)S\A(x)\ 2 becomes dom- 
1 (dt0(x)/4) 2 which only gives negligible corrections to the dispersion of the collective mode 

rf 



inant with respect to n , 

Moreover, in the low-density limit, the coefficient — jr^mp reaches the unity on increasing the interaction 
behaviour and the fact that the mean field value of A 2 , is po/2@, allows us to recognize in the effective action (3.4) 



This 



the contribution — \A(x)\ 2 d t 0(x) where |A(a;)| 2 = po/2 + S\A(x)\ 2 . In these limits (strong coupling and low-density) 
the role of the conjugated momentum to the phase is taken by the square of the modulus of the order parameter: 
indeed, |A(:r)| 2 coincides with p(x)/2. Using the same reasoning of the continuum limit we can then get, as a unique 
equation of motion, a non-linear Schrodinger equation for the field ^ ~ \f\ e%B ■ Since, in this case, 'J coincides with 
the order parameter |A|e lS , the TDNLS results into the Gross-Pitaevskii equation. 

In the weak-coupling limit, the delocalized structure of the Cooper pairs does not permit such a description in terms 
of a Gross-Pitaevskii equation as in a bosonic system. Indeed, we need different equations to govern the dynamics 
of the density and the modulus of the order parameter. On the contrary, when the pairing interaction is strong 
enough to form tightly bound Cooper pairs and the density is low enough that there is no overlapping between them, 
the difference between the square of the modulus of the fermionic order parameter and half of the density becomes 
negligible and they can be identified: |A| 2 = p/2. 



IV. THE PHASE-ONLY ACTION IN THE PRESENCE OF A LONG-RANGE INTERACTION 



We devote this last section to the introduction of a long-range Coulomb interaction Vc(q) = 47re 2 /q 2 in the system. 
In the discussion of the hydrodynamic action in the continuum case we have already implicitly considered a charged 
system when we introduced the electromagnetic field. This time we pay explicit attention to the effects of the Coulomb 
interaction on the structure of the resulting effective action when expressed in terms of the phase only. 

It has been proposed on several occasionally that a low value of the superfluid density and consequently of the 
phase stiffness, enhances the role of the phase fluctuations. It is therefore relevant to study the interplay between the 
phase fluctuations and the low- lying energy modes once the amplitude is kept constant. 

To approach this problem let us consider the electron action (3.1) in the presence of the Coulomb interaction: 
~1 Sir a' Si j Vcifi ~ ^j)4>i ( x i) < t > a(xi)4>t' i x j) < Pcr' ( x j)- We decouple the long-range interaction in the Hartree channel 
and then follow the same procedure as before. Finally we obtain the effective action: 

S eff (p,\A\,6) =- J dtJ2p(%i)dtO(xi)/2 + J dtY,(Eo m )/8(l^co S (9 l -e j )) + 



(ij) 



n P p(<z) + \v c {q) 



5p^ q + S\A\ q U AA {q)S\A\^ + 2S Pq U pA (q)S\A\ 



(4.1) 



Notice that the coefficient of the quadratic part in the density fluctuations has been modified by the presence of 
the Coulomb interaction according to RPA prescription. As discussed in Section III, to compare with Ref.(19) and 
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Ref.(33) we consider here the discrete form — Yl{ij) (1 ~ cos(6>i — Oj)) of the gradient term. We recall again that is 
only in the limit t/\U\ <C 1 that we can consistently derive this term. In particular, to get a cosine term which is local 
in time, we need that the characteristic time over which 6i(t) v aries is much larger than the inverse quasiparticle gap. 



By integrating out all the fields but the phase in the action (hi), we get the phase-only action: 



where 



s («> = \J^E i + v c m( q ) " 2e « e -« + J dt ^ {E ° m)/8 (1 - cos( ^ - 6j)) (42) 



2 n pp n AA - u 2 pA ■ 

H(q) reduces to the "short-range" compressibility x m the limit q — > . 

In agreement with Ref.(33), the action is made of two contributions: the quantistic one (oc aj 2 9 q 0- q ), whose 
coefficient is proportional, in the static limit, to the full compressibility in the presence of the Coulomb interaction 
(and it is therefore vanishing as q 2 ), and the classic term (oc cos(8i — 6j)) whose coefficient is the stiffness related to 
the superfluid density (p s /m*). 

The above phase-only effective action is obtained for a model with long-range interaction and local pairing interaction 
without any a priori assumption. At Vc = 0, the action is equivalent to a x — y quantum model in d+ 1 dimensions. 
On increasing the pairing interaction the stiffness decreases due to the enhancement of the effective mass, while the 
compressibility increases since we are reaching the bosonic limit, where the chemical potential does not depend on 
the densitye3. 

For Vc 7^ 0, the low-momentum and frequency behaviour of the quantistic term is completely dominated by the 
divergence of the long-range interaction and the vanishing of the compressibility In this regime the action becomes: 

^E^ 2 32^%) fl (-«)+ / ^E^)/8(l-co S (^-^)) (4.3) 

1 (ij) 

In this case, the sound mode corresponding to Anderson-Bogolubovil@ mode becomes massive with the plasmon 
frequency given by lo 2 — Ane 2 po / m* in the continuum and by lo 2 — —2ire 2 (Eq" l ) in the lattice model. 

We can extend the analysis to a superconducting system which is coupled to a normal metal through the Coulomb 
interaction. The density fluctuations of the metal will screen the Coulomb interaction in the Hartree channel and 
affect the coefficient of the quantistic part. In this case, the coefficient of quantistic part of the effective phase-only 
action is: 



n(?)/ 



(4.4) 



where e(u>,q) = 1 + Vc(q)TL m (q) and H m (q) is the polarization bubble in the metal. This amounts to replace the 
Coulomb interaction (in the superco nduc tor) with the effective RPA Coulomb interaction screened by the metal. The 



low-energy and momentum limit of (4.4) is j ^e( uj, q)q 2 . In the "static" region (oj -C \q\) the dielectric function is 



proportional to 1/q 2 , the quantistic coefficient (4.4) is a constant and does not modify the behaviour of the short-range 

(Vc = 0) action (4.2). On the contrary, in the dynamic limit the dielectric constant is given by e ==> 1 + iiri^^p- and 
the quantistic coefficient depends on the ui behaviour of the conductivity a(u>). 

Assuming that an action like this also holds ^at finite temperature (with the same expressions for the coefficients) 
one recovers the Emery and Kivelson proposaO for a phase-only action describing HTCS as "bad metals" . In this 
case the onset of the locking of the coherent superconductive phase would be affected by the w-behaviour of the 
conductivity. It is tempting to suggest that, near T c , the role of the external screening system can also be played by 
the "normal" component of the superconducting system at finite temperature. 



V. CONCLUDING REMARKS 



The procedure we have followed to derive the effective (hydrodynamic) action for a superconducting system at 
T = strongly relies on the introduction of the physical fields p, A, A + beside the auxiliary fields p HS , A HS , (A HS ) + 
introduced to decouple the interaction by means of the Hubbard-Stratonovich transformation. The resulting effective 
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action is a proper starting point for selecting relevant features in systems that deviates from the standard BCS 
superconductors. 

The introduction from the outset of the physical density makes it transparent in the action that p{x) plays the 
role of the field conjugated momentum to the phase. As a direct consequence the equation of motion for the field 
*0) = v / p(x)/2e w W has the form of a non-linear Schrodinger equation. However, besides this equation, in general, 
also the equation of motion for the order parameter amplitude must be considered. 

Our procedure allows to follow the evolution of the density fluctuations as the pairing interaction is increased, 
towards the strong-coupling limit when the Cooper pairs can be considered as true composite bosons. In this case the 
density fluctuations are identified with the square amplitude fluctuations and in the low-density limit the two fields 
p(x)/2 and |A(ir)| 2 coincide. As a consequence there is no need of two different equations of motion for p(x) and 
| A (a;) | as in the BCS case, and the low-energy behaviour can be described by a Gross-Pitaevskii equation for the field 
ty(x) which in this case coincides with the order parameter A(x). 

The inclusion of the long-range Coulomb interaction is straightforward and the modifications of the action follows 
the RPA prescription. Integrating out all the fields but the phase we get the_phase-only action in agreement with 
Ref.(33). In this case a comparison with the proposal of Emery and Kivelsontj indicates that their effective phase- 
only action with the coefficient of the quantistic term proportional to the normal conductivity can be obtained if the 
superconductor is coupled to a metallic component. 

This work has been supported by INFM "Progetto di ricerca avanzata" 1996. 



APPENDIX A: THE EFFECTIVE ACTION 



Here we report the derivation of the effective action (1.3 ) in terms of the physical fields p(x),\A(x)\ and 9(x), 
following the steps anticipated in the text. 

First, after having decoupled the interaction with the HS fields p HS and A HS we introduce the physical fields p and 
A in the action by two Gaussian transformations of the quadratic terms (p HS ) 2 and (A HS ) 2 appearing in (1.2). In 
the Nambu formalism, the resulting fermionic-bosonic action, containing both the HS and the physical fields, reads: 



S(r], r?+ , p HS , p, A HS , (A ffS ) + A, A 4 



dx < rj + {x) [idtcra] rj(x) + r/ + (x) 



2m 



+ fx - p HS (x) 



v^ix) + p(x)p HS (x)- 



-r] + (x)7ieA HS (x)air](x) + rj + (x)lmA HS (x)a 2 r](x)+ 



-A(x) (A HS (x)) + + A+{x)A HS (x) - ^P 2 (x) - U\A(x)\ 2 



(Al) 



where r/ + = is the Nambu spinor, <7q is the unit matrix and <Ji is the Pauli matrix vector. Notice that, 

integrating out p HS (x) we get the constraint S(p(x) — rj + (x)a-zr)(x)) that specifies the field p(x) as the local density. 
Analogously, acting on A HS , we obtain 6 (iZeA ~ ^(<f>±<f>i + <f>^(f)^\ and 6 (imA + |(</>;</>| — (j)*<p~^)j . 

Next, we get the action directly in terms of the modulus of the physical order parameter and its phas e. T o this aim, 
after writing A = | A|e ie ( x ) we make the gauge unitary transformation 4>{x) — * <p{x)e l6 ^/ 2 in the action ( Al). Then, wc 
remove the dependencies from the phase in the off-diagonal terms and in A(x) HS \ A(x)\e~ t9 ^ + \ A(x)\e t9 ^ (A HS (x)) + 
by defining the new auxiliary field: A HS (x) — A HS (x)e~ l8 ( x \ Such procedure is convenient in the low-energy limit, 
since it explicits the dependence on the gradients of 6(x). 

After these manipulations the action reads: 



S(r),r) + ,p HS ,p, \A\,6,A HS ,(A HS )+) = / dx* < r] + (x) 



2 m 



a r)(x) 



+r) + (x) 



2m 



p - (d t 6(x)/2 + (x/6) 2 /8m) - p HS {x) 



a 3 r)(x) + p(x)p HS {x) 
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Xj 2 

+ir] + (x)— — <JoT](x) + i] + (x)XmA HS (x)a2r](x)+ 
4m 

- V +{x)TZeA HS (x)cr 1 r 1 (x)+2\A(x)\TZeA HS (x) - ^p 2 (x) ~ U\A(x)\ 2 

By integrating out the fermionic fields we obtain the following bosonic action: 

S(9, p HS , p,\A\, A HS ,(A HS )+) = -iTrlnG- 1 
U 



(A2) 



dx q 



p(x)p HS (x) - —p 2 (x)+2\A(x)\KeA HS (x) - U\A{x)\' 



(A3) 



Here G is the matrix operator representing the one-particle propagator in the presence of the external bosonic fields. 



G~ 



2m 



2m 

.w 2 o 



+ M - P l 



dtO/2 



M) 2 



8m 



03+ 



+ i 



4m 



(T - fteA HS (»o-i - lmA HS (x)a 2 



(A4) 



We re- write G 1 = G 1 (1 — GqT,) by introducing the saddle-point one-particle propagator Go and the self-energy 
matrix S. One has : 



Go" 1 - id t a Q + + fJ> — pg s )* 3 - Ag s a 1 , 



1 2m 



which, in the Fourier space, leads to 



Go — 



Go F 



w + €k A 

us 2 - El V A w - e/c 



with Go(uj 7 k)—G (—u>, —k) and F (u, k)=F$ (uj, k) . Here Aq S = £/|A |, p^ s = Upo/2 are the parameters that 
extremize the action (|A3j). The "self-energy" matrix is given by: 



= <^ 6p HS (x) + 



d t O{x) 



(VO) 



2 1 



8m 



^3 



4m 



2m 



<5^eA HS (x)cti - SXmA HS (x)a 2 



(A5) 



We take X as an expansion parameter, by assuming, as usual, small fluctuations of the fields around their saddle point 
values and small variations of the phase in spa ce and time, since the dependency of the action in 6 is only through 
its gradients. We then rewrite Tr InG -1 in @) as Tr In Gq 1 - Tr£„i(GoX) n where the factor (G S) n is a product 
of propagators connected into a closed loop by the trace, with n insertions of the fluctuating fields. 

To derive the low-energy effective action we retain in ( |A3| ) the first and the second order of the expansion in (Gq£) 
around the superconducting saddle-point: 



S eff (6p,6p HS ,6\A\,6,6KeA HS ) = - / dx*p 



d t 9{x)/2 



8m 



+ 



<lq A r u 



(2tt) 4 [ 4 



-Sp g S P - g + 5p H J q 8p q - U8\A\ q 5\A\_ q + a|A|,afteA« 



[^(?)-D(9)] 



« + l 2' 
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-2i 



1 f Al 



2 9 8m J (2tt) 
d 4 qi 



2 9 8m J (2tt) 4 



4 9l(9l +<l)9q 1 6-qi -q 
31(91 + q)OqiO-qi-< 



[B q + B- q ] SKeA™ 
[B q - B_J 51mA H J„ 



-iSlmA" s [C q - C- q ] SUeA™ - SImA" s 



D. 



■1 ~ + ^~ q ' 



SZmA™ 



+ 5KeA« s 



Dq + -{C q +C^q) 



6KeA HS r 



where the coefficicntsEa are: 



A(q) 



dk 4 



[G (k + g)Go(fc) + G (fc + g)Go (fc)] 



2 J (2tt) 4 

1 /• dfc 3 e+e_ - D+D_ + D_ 

^27 (2^)3 £^EL (D+ + D-) 2 - w 2 



7 /" dk 4 

B (<l) = n I 7^7 [Mk + q)G Q (k) + G "(fc + g)Fo(fc)] 



2 7 (2tt) 4 

Ao f dk 3 
' 4 ./ (2tt)3 



E++E- 



/dk 4 
Go (k + q)G (k) 



dk 3 
(2^)3 



E + E- + e + e_ 



E+E- 

D(q) = l 



(27T) 



(£++£_) 2 



\ 2 



dk 4 

dk 3 1 



F (fc + g)F (fc) 



E+ + E- 



2 7 (2tt) 3 (E+ + E.) 2 - uj 2 



and, in our notation: 



IA 



HS |2 



1/2 



(A6) 



fit = [4 

e+ = (fc + q) 2 /2m — fx — p$ s and e_ = k 2 /2m - fx — p^ s 

Notice that in the q — > limit one gets Go + -Do = —jj — (A) — -Do)- In the BCS limit one recovers: A? — D q 

= N and Bq = 0; while in the strong coupling limit one obtains A — D a ~ — and £? = —(1 — po)Ao^j. 
Since we are interested in the low-energy regime of the system in the derivation of 



have neglected contributions that gives corrections that go beyond this regime. We have indeed ne- 



glected terms associated to the expansion of 



4m 



do whose contribution is proportional 
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-3-0 

2m 17 -?' 



The coefficient of this term 



to J d 4 q^9 q J d 4 k [(fc + q)G (k + q)kGo(k) + (k + q)F Q (k + q)kF (k) 

(J d A k (k + q)Go(k + q)kGo(k) + (k + q)Fo(k + q)kFo(k) ) represents the paramagnetic current and it vanishes in 

the static limit at T — 0, while, at finite temperature, it gives the contribution that, due to the presence of the normal 
component of the superfluid, reduces the stiffness. At zero temperature the first contribution to the action coming 
from this term is of order (\/ 2 9) 2 and then negligible in the low-energy limit. This also shows that the expansion of 

~ + i V a o does not generate time derivatives of 9. This expansion can only give time derivatives of 

(V 2 #), the first non vanishing term being (y 0) . 

The last step consists in integrating out the HS fields to get the effective action ( |l.3[) . Since only the first and the 
second order terms in the expansion in powers of (GqT,) were retained, the action ( |A3[ ) is Gaussian in the auxiliary 
fields even though it contains higher order corrections in the phase fields. We can therefore straightforwardly integrate 
the HS fields to get the action in terms of the physical fields reported in (1.3). 

Notice that in this integration we can ignore contributions coming from the coupling of SX mA" s with the other 
HS fields. The corrections coming from this coupling to the coefficients of the effective action (1.3) represent higher 
order corrections to the hydrodynamic limit. 

We devot e th e remaining few comments on Galilean invariance. Notice that by making explicit the dependence of 
the action (A2) on the phase of the order parameter, the Galilean invariance of the action becomes transparent. A 
Galilean transformation is described by the following change of variables: 



t' 



vt 



v'= V 



where v is the velocity of the moving frame. Under a Galilean transformation A(x) A(x)e 2t \- mv ' r ™v 2 t/2] an( j 
therefore the phase changes into: 

9' \x) — 9{x) — 2mv ■ r + mv 2 t 



In the action (A2) we can individuate two Galilean invariant terms: 

i) idt + *2m V i that can be interpreted as the "covariant" derivative, 

ii) [dt9(x)/2 + (v#) 2 /8m], where dt9(x)/2 is the local chemical potential shift while (\/9) 2 /8m will give rise to the 
superfluid contribution to the kinetic energy of the system, once we define the superfluid velocity v s = (xj9)/2m. 

Notice that, deriving the effective action flA6|) we have spoiled the "covariant" derivative idt +i?^V contained 



in the one-particle "propagator" (A4). As a consequence, for a given order in the perturbation theory, the Galilean 
e effective action can only be recovered taking into account appropriate terms from higher order 



invariance fo 
contribut 
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